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Abstract 

The properties of heavy quark systems change if they are placed 
in a medium other than the low energy vacuum. In a hot Quark 
Gluon Plasma J/\t particles will melt and not exist as resonant states. 
T's, however, because of their smaller size and the dominance of the 
Coulomb potential, will still form as QQ bound states but their prop- 
erties will shift. In particular the T mass may be shifted upward 
by over 150 MeV. If such excited states manage to escape from the 
plasma as a superheated particle, they may serve as a diagnostic of 
the plasma in which they originated. We propose that such a scenario 
is possible and that hot T's will form and escape, thereby providing 
us with crucial information about the Quark Gluon Plasma. 
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1 Introduction 



As the hunt for the Quark Gluon Plasma (QGP) heats up, the search for a 
clean diagnostic signal becomes pressing. The properties of heavy Quarko- 
nium have figured prominently, and promisingly, enal for the formation 

of the QGP. Over a decade ago the melting of the J/^> particle in a plasma 
was suggested M cLS cL si gnal for QGP. However, the subsequent observation 
of J/\I/ suppression in heavy ion collisions was the result of nuclear absorp- 
tion rather than melting [for a review of the field, see 0]. More recent 
observations 0are indicative of a melting effect j|, §, || although, again, the 
interpretation is not unambiguous J7L ^|, 0|. 

The advent of RHIC and LHC may give us the opportunity to clearly see 
this melting effect. The higher energies available to collide larger nuclei will 
produce hotter, more dense material. In addition we expect production of 
heavy Upsilons [see, e.g., [pLOfl l - 111 this note we propose that the properties 
of these heavy Quarkonium may expose the existence of the QGP and help 
in exploring its properties. 

The original suggestion for quarkonium melting in QGP was based on 
the idea that the confining inter quark potential would be screened in the 
plasma and that sufficient screening would liberate the quarks. While re- 
cent discussions of J/ty suppression rely on a more microscope description 
of gluon quarkonium interactions, the original picture is more immediately 
transparent and suggestive of further effects. Our discussion will be framed 
in the context of a macroscopic picture of the QGP, and its inherent screening 
properties. 

Fragile quarkonium such as J/^f, are readily disrupted by screening. The 
more robustly bound T will maintain its integrity as a bound state in the 
plasma but will have its properties, such as its binding energy and hence 
its mass, shifted. Since the screened T will still be small, comparable in 
size to the unscreened J/^f, it will not be highly interactive with the hot 
hadronic gas that will eventually replace the cooling plasma. It thus has a 
reasonable chance of escaping relatively unscathed from its fiery surroundings 
and decaying in a detector with a mass different from its unscreened, normal, 
brethren. This mass shift will be a noticeable signal, and possible diagnostic 
of the Quark Gluon Plasma. The remainder of this note will be a justification 
of the above scenario. 
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2 Hot Upsilon 



2.1 The Screened Cornell Potential 

Debye screening in a conventional plasma is well established. Analogous 
effects are expected in the QGP. A heavy QQ in the medium of the QGP will 
still experience an attractive potential but both the short distance, Coloumb, 
and the long distance, linear, pieces of the potential will be screened. If 
we assume the zero temperature, zero baryon density, interquark potential 
is described by the phenomenological addition of the long range confining 
potential and short range Coulomb potential we have the Cornell potential 

V = ar-- (1) 

r 

We expect that the Coulomb, perturbative piece will undergo Debye 
screenning and behave like 

- -e"", (2) 
r 

with /i(T), the screening mass, an increasing function of the temperature T. 
The Debye screening radius is . The linear part of the potential should 



also be screened and we adopt the suggestion of [11 



-(1-e-n, (3) 
leading to the screened QGP potential 

V(r, u) = -(l- e-n - -e-» r (4) 
fi r 

Evidently the confining character of the potential is gone, just as we 
anticipated. However, we still have an attractive potential which can bind 
a QQ pair into quarkonium if the quarks are sufficiently massive. It is a 
dynamical question which states will bind and which will not. Even for 
moderately high values of T, the screened potential (4) will bind bb into a 
bound "hot" T with properties different from the properties, determined by 
(1), of the T at T = 0. This contrasts with the case of J/^f which readily 
loses its binding and melts. This was the initial motivation for the proposal 
that J/ty suppression would be a signal for the QGP. The fact that a modified 
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T may survive offers the opportunity to study some properties of QGP in 
more detail. 

In order to study the modifications introduced by screening we need an 
acceptable T = potential. By now this is a routine task and has been done 
numerous times. We use here the values presented in [[12] (based on a fit to 
the spin-averaged energy levels of T family): 

a = 0.470, a = 0.186Gev 2 , m b = 4.753Gev (5) 

We now ask what happens to the bound states when the screening turns 
on as T increases. In our simple model all the T dependence is contained in 

2.2 The T — dependence of // 

The T-dependence of the Debye mass is a subject of intense ongoing 
theoretical research [13-18]. A survey of the literature leads to the conclusion 



that the leading-order perturbative result |L3 



V LO (T)/T = (N c /3 + Nf/Q) 1 ' 2 ■ g(T) (6) 

consistently underestimates the lattice results for /i(T)/T. Significant non- 
perturbative contributions are present up to very large temperatures. (In 
eq.(6) N c is the number of colors, Nf is the number of light quark flavours 
involved and g(T) is the T-running coupling constant.) 

We are interested in the behaviour of the Debye screening mass at tem- 
peratures ranging from one T c to a few (4-5) T c {Tq is the deconfinement 
transition temperature). In this range of temperatures there are lattice re- 
sults for both SU(2) [|14], [HJ and SU(3) pure gluonic theories. 

In the case of SU(2) both |14j and Jl5| report a practically linear In- 
dependence of /i, albeit with different slopes: fi(T) ~ 2 • T from [jrjj and 



is 



/i(T) ~ 2.8 • T from [|T^]. For the SU(3) case, the result reported in [[16 
/i(T) ~ 6.3 ■ T. 

Results including dynamical quarks are not yet available from lattice 
measurements. Nevertheless, model calculations based on the method of 



dimensional reduction of high T QCD (see for instance [|l7|j, ||18|| ) offer some 
guidance for the dependence of the screening mass on the number of quark 
flavors. Ref. |T7| calculates the value of the ( electric) screening mass at 
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T = 2T C , for SU(3) and three flavors, to be n(T) ~ 6.6 • T. (The same work 
offers the screening mass for pure SU(3) at T = 2T C namely /x(T) ~ 4.6 • T). 

In principle the lattice measurements should be the most reliable ones 
as they use a fully non-perturbative method. Therefore we use here the 
previously mentioned linear T-dependence of /i allowing a range for the pro- 
portionality constant: /x(T) ~ (2.5 — 6) • T. As a working case we choose 

Ai(T) = 4 • T (7) 
2.3 The Hot Upsilon 

We are most interested in the mass and the size of the possible hot bound 
states. The mass shift will serve as a marker for the unusual environment 
in which the bound state forms while the size will be a useful guide to the 
ability of the shifted, superheated state to slip through the hot hadronic 
environment and escape. 

We solved the Schrodinger equation numerically with the potential given 
by eqs.(4),(5). Table 1 shows the masses and sizes of T(IS') states for several 
values of the screening mass \i. 



/i(Gev) 


Mass(Gev) 


radius (fm) 


0.00 


9.454 


0.23 


0.36 


9.555 


0.26 


0.72 


9.615 


0.31 


1.08 


9.639 


0.46 


1.40 


9.635 


1.19 



Table 1: Mass and radius of T(IS') for different screening masses. 

The mass of the T(IS') is seen to increase as \x (and T) increases from 
to lGeV. The size of the state also increases but remains relatively small. 

For instance, at \i = 0.72GeV, corresponding to T = T c = 180MeV, 
the hot T(15) has a mass of 9.615GeV (159MeV^ greater than the T = 
mass) and a size of 0.31/m, significantly smaller than that of a cold J/^> 
(~0.45/m, 0). 

Beyond \i ~ 1.50GeV the potential of eq.(4) does not allow bound states, 
i.e. the states "melt" . 
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The mass shift is easy to understand. The bb quarks in the T primarily 
feel the influence of the Coulomb part of the potential. They are so close 
together that they are hardly aware of the confining part and thus barely 
notice its disappearance at high T. It is this insensitivity to the confining 
potential that is responsible for the survival of the T at high T. If the 
Coulomb potential is screened, its negative contribution to the bound state 
energy is weakened, increasing the T's mass and radius. We note here that 
these (relatively) small changes in the bound state properties will not affect 
significantly the decay width. The lifetime of the hot state exceeds, by large 
factors, any strong interaction timescale. 

If this superheated state escapes and can be observed it will indicate the 
formation of a macroscopic medium with a screened potential. The most 
optimistic scenario would allow this shifted mass to serve as a thermometer 
of the original QGP. Of course, at very high T even T(1S) will melt, but, as 
we shall argue, in the temperature range of RHIC we expect mass-shifted Ts 
to form. 



3 Fate of Hot Upsilons in Heavy Ion Colli- 
sions 

3.1 Hot States and P T Dependence 

We begin by looking at a heavy ion collision which produces a QQ system 
at t — 0, at rest in the lab (CM) frame. After a time t , a large (effectively 
infinite) QGP, at temperature T , forms. The QGP will begin to cool. We 
assume the cooling rate is related to the expansion time as given by longitu- 
dinal isentropic expansion equation [I9| 



s(t )*o = s{t)t (8) 

where s(t) denotes the entropy density at time t. 

There are three temperatures that are of significance: 

a) To is the initial temperature of the plasma. 

b) T c is the transition temperature above which a QGP can form. Recent 
results from lattice computations suggest T c ~ 170 — 190MeV |2(J. In this 
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paper we will use the value T c = 180MeV\ Below T c we have a hot hadron 
gas. 

c) We define T m as the melting temperature above which bound QQ 
states can no longer form. In our model it is the temperature for which the 
potential (4) no longer supports bound states. For the T(IS') the melting 
occurs for a screening mass /i m ~ 1.50GeV. If we use eq.(7) this corresponds 
to T m ~ 375MeV. 

The corresponding times are also of significance: 

a) to is the time of formation of the plasma. 

b) t c is the time at which the plasma cools below T c and becomes a hot 
hadronic gas. 

c) t m is the time at which the plasma cools below the melting temperature 

T 

These times acquire their significance when compared to the formation 
time tf, the time it takes for the initial QQ pair, at rest, to form a quarkonium 
bound state. A simple model for tf [22] is the time required for the two 



quarks, forming a specific quarkonium state, to separate a distance R equal 
to the size of the corresponding state. The velocity of separation v is taken 
to be the average radial velocity of the quark as given by the potential model 
for that state, so tf = R/v. 

The formation time tf increases with the temperature T, since R increases 
and v decreases when T increases ( the screening mass \x increases and the 
inter-quark potential weakens). We will be interested in what follows in tfo 
and tf c , the formation times at T = and T = T c , respectively (note that 
tfo < tf c ). 

We are now ready to discuss the fate of heavy QQ bound states. 
A. For T c < T m . 

1. t c < tfo m . the quark- ant iquark pair escapes from the plasma and forms 
a normal, T = 0, cold state within the hadronic gas. 

2. tfo < t c < tf c : the pair does not reach the minimum size necessary 
to form a hot bound state and the quarks are too far apart to form the cold 
state in question when the QGP is replaced by the hadronic gas. The state 
melts, although under some limited circumstance an excited, cold resonance 
can form (for instance, T(2S') instead of T(1S)). 
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3. tf c < t c : the pair forms a bound state in the plasma, a hot state. Its 
properties (e.g. mass and wave function) will be determined by the screened 
potential appropriate for some T, with T m > T > T c . 

B. For T c > T m the picture is simpler: 

1. t c < tfo'. the pair forms clS db cold quarkonium state; 

2. t c > tf Q : the quarkonium state cannot form in the plasma (too hot) 
and it melts or may form a cold, higher resonance. 

Thus, for QQ at rest, the formation of hot quarkonium states is deter- 
mined by the conditions: 

t fc <t c and T c < T m (9) 

Most QQ are, however, not produced at rest with respect to the QGP rest 
frame. They will have transverse momentum P T and we must account for the 
relativistic time dilation when computing the relevant tf. The formation time 
in the QGP rest frame for a system of mass M and moving with momentum 
Pt is 

<'/(Pt)=^1 + ^ (10) 
Hence eq.(9) is replaced by 

t' /c (P T ) < t c , T c < T m (11) 

In order to relate the times t m and t c to temperatures we need an ex- 
plicit expression for the entropies appearing in Eq.(8). As usual we take the 
simplest reasonable form, i.e. the entropy for a free gas: 

s(T) = const ■ T 3 (12) 

Combining Eqs.(8),(9),(10),(ll) and (12), we find that those QQ pairs 
produced with transverse momentum in the range 

(£) < (£) (|) - 1 < l3 » 

will form hot quarkonium states. 
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3.2 Hot T States at RHIC 

Let us use eq.(13) for T(IS'). This state melts at \i m ~ X.^GeV and 
the corresponding melting temperature is T m = 375MeV, according to our 
choice /i(T) = 4 • T. The formation times needed are t/o = 0.76 fm and 
tf c = 1.32/m. 

We will use values for t and T based on the expectations at RHIC. 



Recent estimates |21j suggest the QGP plasma should form at to — 0.17 fm 
with an initial energy density eo — 98.2GeV/ fm 3 and initial temperature 
T = Q20MeV. 

Using these values together with M = 9.615Ge^ in eq.(13), we find that 
T(IS') states with in the range 

< P T < AVGev (14) 

will be produced hot and have a mass shift of about 160MeV. Thus, accord- 
ing to our model, essentially all T(IS') produced at RHIC will be "hot"! The 
question we now ask is, will they survive? 

3.3 Survival 

If hot T's form in a QGP, as outlined above, they will be interesting only 
if we know of their existence. The most promising signal of this existence is 
if some of the Upsilons escape from the plasma and its surrounding fireball. 
Since the hot T is still small, a geometric picture of scattering would imply 
a reasonable likelihood of escape. 

We picture a large, spherical plasma wich cools down to a hot hadron gas 
at time t c . At this time all the hot Ts will be immersed in the hadron gas 
where they can scatter and cool down or be disrupted. We will quantify this 
in a simple absorption model, following the treatment used by Karsch and 
Satz f22j for the absorbt ion scenario in an expanding, hot hadronic gas. 



The survival probability for the T entering the hadron gas at time t c and 
surviving until time r is 



S(t) = exp — / dtn(i)o (15) 

L Jtr. J 



Here n(t) is the density of the medium, a is the appropriate hadronic cross 
section for Y(IS') and r is the time it takes for the hot hadron gas to become 
so cool and dilute that no further interactions are likely. We take this 'freeze 



9 



out' conditions to occur when the energy density is e(r) ~ 0.3GeV/ fm 3 . 
If we continue to assume isentropic longitudinal expansion for the hadron 
gas (treated as an ideal gas) as well, the relation between densities at the 
relevant times is 

n(r)r = n(t c )t c (= n(t )t ) (16) 
and the survival probability becomes 

S(t) = — — , k = n c t c a = n t a (17) 
n(t c ) 

As a final step, we use the relation between the density n and the energy 
density e for an ideal gas, n = (|e) 3//4 . This allows us to express S in terms 
of the "accessible" quantities eo,io an d cr: 



S 



VH fee (18) 



Using T c = 180MeV in the energy density for an ideal gas of gluons and 
3 flavours of massless quarks 

we find e c ~ 2.1AGeV / fm 3 . Then, using the above-mentioned value 
e(r) — 0.3GeV/ fm 3 , the survival probability is: 

S(r) = [o.i4]( 3 / 4 )-( 2 / 3 ) 3/4 - e ° 3/4 - to - (19) 

The crucial unknown in this equation is the cross section for a hot T 
interacting with hadrons.The escaping hot T(IS') can suffer two types of 
absorbtion in passing through the hot hadronic matter: 

a) It can be dissociated by hadronic collisions, leading to supression of the 
T signal similar to that initially observed for J/^f. The survival probability 
is given by eq.(20) with a the inelastic cross section. 

In the limit of large quark mass, QCD allows us to calcculate [23-26] the 
inelastic scattering cross sections of J/ty and T on light hadrons (pions or 
nucleons). The asymptotic cross sections are estimated to be = (2— 3)mb 
and (Tx-tt = (0.6— 0.8)mb, which are comparable to the geometric cross section. 
However, the asymptotic cross section is approached very slowly because of 
the presence of a strongly damped threshold effect [23-25]. For instance, 
for Er < lOOGeV, ox^ < O.lmfe. The corresponding predictions for J/ty 
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inelastic scattering have been succesfully tested against photo-production 
cross sections [27]. 

Therefore we adopt a (conservative) estimate of 0.3m6 for the inelastic 
(dissociation) cross section of the hot T(IS'). 

Using this value in eq.(19) in conjunction with cq = 98. 2G eV / fm 3 and 



t = 0.17/m PU, we find S ~ 0.56 i.e. 56% of the produced T(1S) still 



survive the hadronic fireball. 

b) T(15) can scatter quasi-elastically into the normal, T = 0, T(1S). 
This is the process that is of most interest to us, since it decides wheather 
or not a significant number of hot T(15) could be detected. 

To make an estimate to the fraction of hot T(1S) versus cold (T = 0) 
T(IS') we therefore need to know the cross section a qe associated to the 
reaction 

T(lS)(hot) + hadron -> T(1S) (cold) + hadrori (20) 

The QCD approach to heavy quarkonium scattering can also be applied 
to the elastic and quasi-elastic scattering |23|, ^6[. These cross sections are 



expected to be quite small, considerably smaller than the inelastic cross 
sections, an expectation that is again born out by analysis of the photo- 
production data [27]. 

To be conservative we use a value a = O.lmb in eq.(19) and obtain S — 
0.82, i.e. 82% of produced hot T(1S) states will survive their traversal of the 
hot hadron g clS clS hot, mass shifted states. 

To conclude this section, let us review the results: 

i) all directly produced Y(IS') states with transverse momenta below 
A9GeV are hot, with a mass shift of over 150Me^; 

ii) 56% of them will survive the passage through the hadronic aftermath 
either as hot or cold states; 

iii) the fraction of surviving T(1S) which are superheated is (at least!) 
82%. 

3.4 Sensitivity to input parameters 

Our result is sensitive to the values we have assigned for the initial energy 
density eo(= 98.2GeV/ fm 3 ), the hadronization temperature T c = 180MeV^ 
(wich determines e c ), the formation time of the plasma to(= 0.17/m) and 
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the quasi-elastic cross section a(— O.lmb), none of wich are precisely known. 
Our key result, the fraction S of surviving T's that are superheated, changes 
as follows when we vary these values one at a time (within reasonable limits): 

a) e = 80GeV/fm 3 gives S = 0.85 and e = 120GeV/fm 3 gives S = 0.80 

b) T c = UOMeV gives S = 0.84 and T c = 190MeV gives S = 0.81 

c) Larger value of a = 0.3mb gives S = 0.56 

d) varying t affects twofold our result for the fraction of hot Ts available 
for detection. Firstly, to enters eq.(13) which determines the maximum Pt 
for hot Ts (according to eq.(13), to has to be larger than 0.032/m in order 
to have hot T's formed). Secondly, t enters eq.(19), affecting the value of S. 
Values of t from 0.1 fm to 0.25 fm ensure values of values of Pxmax between 
28GeV and 73GeV and values of S between 89% and 75%. 

This range of parameters is probably sufficiently broad that it might 
encompass changes in the underlying models we have used to describe our 
effect. 

3.5 Contribution from Higher States 

Our discussion dealt so far with directly produced T(IS') states. We look 
now into the contribution from the excited states. 

In Table 2 we present the formation times (at T = 0) and melting fi 
values for the higher bottomonium states that contribute significantly to the 
T(IS') production by cascade decays. 





T(2S) 


r(3S) 




X(2P) 


radius at T=0 (fm) 


0.51 


0.75 


0.41 


0.66 


format ion t ime t j ( f m) 


1.81 


2.46 


2.50 


2.88 


melting jj, (GeV) 


0.63 


0.42 


0.58 


0.41 


T m [from eq.(7)] (GeV) 


0.16 


0.11 


0.15 


0.10 



Table 2: Bound state properties of excited T states 



We notice that the melting temperatures, according to eq.(7), are below 
T c and therefore these states cannot form in the QGP. Their formation has 
to take place in the post-plasma era (that is, within a hadronic medium) and 
consequently their binding properties should be the normal ones (there are 
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no hot Xb% e.g.). 

Recalling the analysis of 3.1, we determine the minimum transverse mo- 
mentum these states should have in order to escape melting within the QGP 
at RHIC (T = 620MeV). 





T(2S) 


T(3S) 


Xb(lP) 




PTmin(GeV) 


37 


27.3 


25.7 


22.5 



Table 3: Minimum tranverse momenta to escape melting in QGP (RHIC) 

We find that the feed-down production of T(IS') from higher states con- 
tributes only for transverse momenta larger than the values shown in Table 
3. Consequently, we expect that at RHIC all T(IS') with P T below 22.5GeV 
are directly produced. 

3.5 Discussion on functional dependence of /x(T) 

The numerical results presented so far assume the validity of eq.(7). 
Changes in functional dependence fJ>(T) generate changes in the parame- 
ters tf c (the formation time at T = T c ) and T m (the melting temperature). 
According to the analysis presented in 3.1 we find the following results when 
the functional dependence fi(T) is allowed to vary from 2.5 • T to 6 • T: 

i) For n = 2.5 • T, all T(1S) with transverse momenta < Pt < 66GeV 
are formed as hot states. The only change regarding the higher states is that 
now there are hot T(2S) formed with transverse momenta < P T < 5.8GeV 
that will eventually contribute to (cold) T(IS') signal. 

This possibility will disappear at /i ~ 2.8 • T and for \xjT = 2.8 — 6 there 
is essentially no change, i.e. the conclusions of 3.4 are entirely valid, as far 
as the contribution from the higher states are concerned. 

ii) When /i/T increases to 6, the momentum range of hot T(IS') states 
decreases to < P T < 24.8GeV. 

Note that in the (simple) model we have used to find the survival proba- 
bility of a hot T(IS') state, the function S (see eq.(19)) has no Pt dependence. 
Thus there are no changes in S when /x(T) is varied as mentioned above. 

We then see that despite the large range allowed for the fJ>(T) functional 
form, a significant number of hot T(IS') states could be produced at RHIC 
with negligible feed-down contribution from higher resonances. 
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4 Conclusions 



We conclude that there is a reasonable possibility that hot T(IS') states, with 
masses measurably different from the T = value will be produced, escape 
and be detected. More specifically, our prediction is that about 80% of all 
T(1S) states detected at small rapidity and with transverse momenta below 
25GeV will have a mass shift of ~ 150MeV, this shift being caused by their 
formation in the screening QGP. 

The arguments we have presented are partly qualitative, similar to the 
original arguments in favor of J/^ suppression. Nevertheless the strong 
physical pictures that underlies these arguments leads us to propose a search 
for superheated Ts as a signal for Quark Gluon Plasma formation. 
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